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1 Introduction 

The purpose of this paper is to describe algorithms for computing various functors for 
algebraic P-modules, i.e. systems of linear partial differential equations with polynomial 
coefficients. The algorithms enable us to perform actual computations (with limitation 
caused by the complexity) by using e.g. a program kan ||3^ developed by the second 
author, as well as to establish theoretical computability of some fundamental functors in 
the P-module theory. 

Let K be an algebraically closed field of characteristic zero and let X be the affine 
space K"' with a positive integer n. We denote by Ox and Vx the sheaves on X of rings of 
regular functions and of algebraic linear differential operators respectively (cf. Bernstein 
[|], Bjork @], Borel et al. 0). Let M and M be coherent left Px-modules. 

Various functors are defined for (especially for holonomic) D-modules and play the 
fundamental role (see and also e.g., Kashiwara |]14|, |]16|, Mebkhout ||2^ for their 



analytic counterparts). Among such functors, we are concerned with the following: 

1. The cohomology groups of the restriction Jviy ■= Oy ®Ox of to F as left 
Dy-modules, where y is a non-singular subvariety of X and (8)^ denotes the left 
derived functor (cf. [jl2|) of the tensor product. 



2. The cohomology groups Ext]y^{J^, K[[xi, . . . ,Xn]]) with coefficients in the formal 
power series solutions of Ai , which equal to those with coefficients in the convergent 
power series solutions if Ai is regular holonomic (cf. Kashiwara-Kawai fL^). 

3. The tensor product J^®Ox-^ S'^d, more generally, the torsion groups Tor {M.,J\f), 
as left Px-niodules. 

4. The localization M[f'^] := ®Ox of as a left Px-module, where 
/ G K[xi, . . . , Xn] is an arbitrary non-constant polynomial. 



5. The (algebraic) local cohomology groups 7i[yj(A^) with support Y as left Vx- 
modules, where Y is an arbitrary algebraic set of X. 



It was proved by Kashiwara that these are all holonomic systems (the second one is 
a finite dimensional vector space) if so are Ai and Af. 

Let us remark that ii K = C and Ai is Fuchsian along Y in the sense of Laurent and 
Moteiro-Fernandes p2[, which is the case if Ai is regular holonomic in the sense of 



then there exists an isomorphism 

RHom vx {M, O^") |y ~ RUom VyiMl-, O^"") 

in the derived category of sheaves of C-vector spaces; here O'^ and denotes the 
sheaves of holomorphic functions on X and on Y respectively, and RTiom the right 
derived functor of Horn . Thus roughly speaking, Jviy corresponds to the system of 
partial differential equations which the solutions of M. restricted to Y satisfy. Similarly, 
M. ®Ox corresponds to the system which the product of solutions of AA and of N 
satisfies. 

As was observed by Galligo |lT| and was developed by several authors (e.g. [0, 



3^, l2^, ||2^, [m, [^) the notion of Grobner basis and the Buchberger algorithm p 
are essential in the algorithmic study of P-modules as well as in computational algebraic 
geometry (cf. [|], |^). By using Grobner bases for the Weyl algebra, we give algorithms 
for computing the objects listed above under some conditions on A4. and A/", which are 
certainly satisfied if M. and M are holonomic. These algorithms also apply to the analytic 
counterparts of these functors as long as the input "D-module is defined algebraically. 

We first give an algorithm for the restriction (Algorithm |5.4|) when F is a linear subva- 
riety of arbitrary codimension under the condition that M. is specializable along Y , which 
is the case with an arbitrary holonomic Px-module M.. Here M. is specializable along Y 
by definition if and only if there exists a nonzero 6-function, or the indicial polynomial of 
M. along Y . We also give an algorithm to compute the 6-function (Algorithm |4.6| ). 

Our method consists in computing a free resolution of M. that is adapted to the 
so-called V^-filtration associated with Y . Such a free resolution tensored with T>y^x '■ = 
{T>x)y gives A4y, but it is not a complex of coherent Py-modules in general. Then we use 
information on the integral roots of the 6-function to truncate the complex and obtain a 



complex of finitely generated free Py-modules. The first author gave in |]3T| an algorithm 
for the case where Y is of codimension one without using free resolution. 

This algorithm for the restriction also solves the other problems by virtue of some 



isomorphisms provided by the P- module theory, especially those described in flj]. See 



Algorithm |6.2| for the tensor product. Algorithm |6.4| for the localization, and Algorithm 



77^ for the algebraic local cohomology groups. Finally the computation of the restriction 
for the general case where Y is not necessarily linear reduces to that of local cohomology 
through the so-called Kashiwara equivalence |]14|, which claims the equivalence of the 



category of coherent Py-modules and that of coherent Px-modules supported by Y. 

Algorithms for the local cohomology groups have been given in when Y is of 
codimension one, and by Walther |R3 under the assumption that Ai is saturated with 



respect to Y. An algorithm for the localization has been given in [^] under the condition 



that At is /-saturated. These assumptions are removed in the present paper. As another 
application of the restriction algorithm, we can also obtain an algorithm for integration of 



a module over the Weyl algebra, which enables us to compute the de Rham cohomology 
groups of some algebraic varieties. See for details. 

Finally we discuss how to get the free resolution mentioned above. For that purpose, 
we apply Schreyer's method for free resolution in the polynomial ring (see e.g. 0) to 
the ring of differential operators. In doing so, we need some modification because of 
the non-commutativity and the fact that the term order we use is not a well-order. We 
have two methods to cope with this difficulty: one is the homogenization with respect to 
the V^-filtration by the first author (p8[,|2^); the other is what we call the homogenized 



Weyl algebra which was introduced and implemented by the second author in the 2nd 
version of kan/sml |^8[ that was released in 1994, but has not been published in the 



literature. A similar method was employed by Assi et al. independently and applied 
to the computation of the slopes of a P-module. 

We have implemented the algorithms by using kan/sml for computations of 
Grobner bases and free resolutions in the Weyl algebra, and Risa/Asir for factor- 
ization and primary decomposition in the polynomial ring. 



2 V-filtration and free resolution 

Let K be an algebraically closed field of characteristic zero. We fix positive integers d and 
n. Let X be the affine space with the coordinate system (t, x) = (ti, td,xi, . . . , Xn)- 
We denote by dt = {dt^, ...,dt^) and dx = {dxi, ■ ■ ■ ,dx„) the corresponding derivations 
with = d/dxi, dt^ = d/dtj. We use the notation := a;"^ ■ ■ ■ x"", 9f := d^l ■ ■ - d^", 

:= 9r := dl'-'-d"^: for a = P = (A,...,/?^ G N" and 

/i = (/ii, . . . , fid), ^ = {^1, ^d) £ N'^, where we put N := {0, 1,2,.. .}. We also use the 
notation |a| := ai + . . . + 

Let Y be the ci-codimensional linear subvariety of X given hjY:={{t,x)EX\t = 0}. 
Let Ox and Oy be the sheaves of regular functions on X and on Y respectively. We denote 
by Dx and Vy the sheaves of rings of algebraic linear differential operators on X and on 

Y respectively. 

Let Ai he a coherent left P^-niodule on X. Then the set of the global sections M : = 
T{X,J\4) is a finitely generated left module over the Weyl algebra Aa+n '■= ^{X,Vx)- 
Conversely, for a finitely generated left A^+^-module M, its sheafification 7V1 := Vx^A^+n 
M is a coherent Px-module. More precisely, this correspondence gives an equivalence 
between the category of finitely generated A^+^-modules and that of coherent Px-modules 
(cf. 0,0). Hence we could work only in the first category. However, as to e.g., the 
restriction functor, it would be preferable to work in the latter category since a coherent 
Pjsf-Kiodule can be specializable along some Zariski open subset of Y but not along whole 

Y (cf. Section In any case, actual computations are done for modules over the Weyl 
algebra. 

In the sequel, we define the notion of free resolution adapted to the ^-filtration. Let 
VxIy be the sheaf theoretic restriction of Vx to Y. Let Jy := Oxti + . . . + Oxtd be the 
defining ideal of Y. Then for each integer k we put 

F^iVx) := {PeVx\y\P{Jyy eiJyy-' for any j > fc} 



with the convention Jy = Ox for j < 0. This is called the V-filtration attached to Y (cf. 
I^,p3|). More generally, given an r- vector m := (mi, ...,mr) of integers, we put 



i=l 



where the canonical generators of Vx- We may assume that Ai has a 

presentation Ai = V^/M on X, where A/" is a coherent left "Dx-submodule of P^f- 
fact, M := r(X, A^) can be written in the form M = A'j^^^/N with an integer r and an 
A^+n-submodule of v4JJ_^„. Then M := Vx ^ satisfies the above property. Let Ui 

be the residue class of Cj G Vx in A^. Then for m G Z'", we put 

F^[m]{Af) := AfnF^[m](V'x), 

F^[ni]{M) := F^-"''{Vx)ui + --- + Ft"'^{Vx)Ur 

for each integer A; G Z. The graded ring and modules associated with these filtrations are 
defined by 



gr^(I?x) := ^F^iVx)/Ft\Vx), 
fcez 

gr^[m](P^) := ^ F^[m]iV^x) / F^'l^Wx) 



ryim 



fcez 

W) := ^F^[mm/Ft'[ni]{Af). 
fcez 



Then gry[m](A/') and gry[m](A4) are coherent left gr(Px)-iiiodules. If m is the zero 
vector, we shall omit the notation [m]. 

For a nonzero section P of Vx\y, let k = ordy[m](P) be the minimum G Z such 
that P G F^[m]{V'x). (We put ordy[m](0) := -oo.) Then let ay[m](P) be the residue 
class of P in 



gr^[m](P^) := [m](I)^)/Fr'[m](2^ 



Definition 2.1 Let Ai = P^/A/" be as above. Let us consider a free resolution 

pn-i ^ . . . p^o ^ (2.1) 

of Ai, where ipi are homomorphisms of left Px-modules, and if is defined by V5(ej) = Ui 
for i = l,...,ro with tq = r. This free resolution is said to be adapted (or strict with 
respect) to the V-filtration -Fy[m] if and only if there exist vectors mi G Z^%..., G Z*"' 
such that 

V',+i(F^[m,+i](I?r^)) C F^[m,](I?^) 



holds for j = 0, 1, — 1 with mo = m and that 

F^[mi]m) ^ . . . ^ F^[ni,]{V^^) ^ F^[mo](P^°) ^ F^[m](^) ^ 

is an exact sequence for any k E Z. We call mi,...,m; the shift vectors associated with 
the free resolution 

The definition above is general in the sense that it is local and also applies to the 
analytic case (cf. Section However, from the computational view point, working in the 
Weyl algebra would be more convenient: Let An and A^+n be the Weyl algebras on the 
n variables x and on the d + n variables {t,x) respectively with coefficients in K (cf. [Q). 
Put L := N2('^+") = N'^ X N'' X N" X N". An element P of is written in a finite sum 

P = E E a^^^^rd'^x'^d^^e, (2.2) 

with ttf^^afSi G -ft^, ei := (1, 0, . . . , 0), . . . , := (0, . . . , 0, 1). Put 

F^{Ad+n) := {P= E T.^M^)t''9td^^Ad+n\a^,/3{x)eK[x]}, 

\u\~\n\<k 13 



i=l 

F^[m](M) := + ■ ■ ■ + F^-'-'-(A,+„K 



where Ui is the residue class of in M. 

The following lemma follows immediately from the definition: 

Lemma 2.2 Let p be a point ofY and P a germ ofDx at p. Then P belongs to Fy{T>x)p 
if and only if there exists a{t, x) G K[t, x] such that a{p) ^ and a{t, x)P G FyiAii^n)- 

By using this lemma and the fiatness of Vx over A^+n, we can easily get the following: 

Proposition 2.3 Let 

be a free resolution of M adapted to the Fylm]- filtration, i.e., ipi are homomorphisms of 
left Ad+n-T^odules, ip is defined by f^Ci) = Ui for i = I, ...,^"0 with vq = r, and there exist 
mi G 7/^ xni G 7/'- such that 

z^,+i(F^[m,+i](A^,V'«)) C i^^K](A^,Vn) 
holds for j = 0, 1, — 1 with mo = m and that 

F^V^iW^n) ^ • ■ ■ ^ F^[m,]{AZn) ^ F^[m]{M) 

is an exact sequence for any A; G Z. Under this assumption, the exact sequence ( ^.Sj) 
tensored by Vx from the left gives a free resolution of adapted to the Fylm]- filtration. 



3 Grobner bases and free resolution 



The purpose of this section is to show that Grobner bases homogenized with respect to 
the V-filtration provide a free resolution adapted to the V-filtration. An alternative and 
more efficient method will be described in Section ^. 

We fix a natural number r and a vector m G Z*". Let -< be a well-order (i.e. a linear 
order) onL x {l,...,r} which satisfies 

(a, i) -< {f3,j) implies (a + 7, i) -< (/3 + 7, j) 

for any a, /3, 7 G L and i,j G {1, . . . , r}. (3.1) 

Then we define a total order on L x {1, r} by 

(/i, z/, a, (3,i) -<F (yu', v' , a', j) if and only if 
|z/ — yu| + rrij < |z/' — /i'l + mj or else 

|z/ — yu| + rrij = \u' — /i'l + mj, (/i, i/, a, /5, i) -< (/i', z/', a', j). (3.2) 



Let P be a nonzero element of v4J^_,_„ which is written in the form (|2.2|) . Then the leading 
exponent lexpp(P) G L x {1, . . . , r} of P with respect to -<i? is defined as the maximum 
element of {(/x, i/, a, /5, i) | a^^a/si 7^ 0} in the order -<p. Moreover, for {fi,h',a, = 
lexpp,(P), the leading coefficient of P is defined by Icoef ir(P) := a^yai3i- The set of 
leading exponents Ep{N) of a subset of v4^+i is defined by 

Ef{N) := {lexp^(P) | P G AT \ {0}}. 

Definition 3.1 A finite set G of generators of a left A^+^-submodule A^ of is called 
a Grobner basis of A^ with respect to -<i7' (or an P[m]-Grobner basis), if we have 

Ef{N)= U(lexp(P) + L), 

PeG 

where we write 

{a,i) + L = {{a + f3,i) | /3 G L} 

for a & L and i G {1, . . . , r}. 

We define an order -<h on N x L x {1, r} by 

(A, a, i) -<H (A', a', i) if and only if A < A' or else 

A = A', (A,a,z)^(A',a',j), (3-3) 

where A, A' G N, a;,Q;' G L and z,j G {1, ...,r}. It is easy to see that -<h is a well-order 
and satisfies 

Lemma 3.2 // |z/ — yu| + — A = \v' — + mj — A', then we have 

(A, /i, I/, a, /3, i) -<H (A', /i', z/', a', j) 

z/ anc? only if 



We introduce an indeterminate to which commutes with any element of Ad+n in order 
to define the homogenization. 

Definition 3.3 An element P of y4rf+„[to]^ of the form 

is said to be F[m]-homogeneous of order k if a\^j,pai3i = whenever — fi\ — X + rrii ^ k. 

Definition 3.4 For an element P of of the form (|2.2|) , put 

k := min{|z/ — fi\ + nii \ a^^api 7^ for some 
Then the F[m]- homogenization h{P) G y4d_|_„[to]'' of P is defined by 

Then h{P) is F[m] -homogeneous of order k. 

When m is the zero vector, we simply say F- homogeneous instead of F[m]-homogeneous. 

Lemma 3.5 If P G Ad+„[to] is F -homogeneous and Q G y4d+„[to]'' is F[Yn\-homogeneous, 
then PQ is F[m\-homogeneous. 

Lemma 3.6 For Pi, ... ^P^ G put P = Pi + ■ • • + Pfc. Then there exist /, Zi, . . . , 4 G 

N so that 

t[h{P)=t'^h{Pi)+--- + t'^'h{Pu). 

Let us define w : N x L x {1, . . . , r} — > L x {l,...,r} by ca(A, /x, z/, a, /3, z) = 
(/i, z/, a, /5, i). For a nonzero element P = P(to) of y4rf+„[to]^, let us denote by lexpji^(P) G 
N X L X {1, . . . , r} and lcoef//(P) G K the leading exponent and the leading coefficient 
of P with respect to -<h- 

Lemma 3.7 1. IfPito) G Arf+„[to]'' is F[m]-homogeneous, then we have \expp{P{l)) = 
t^(lexp^(P(to))). 

^. For any P G tfe /iave lexp^(P) = 'a7(lexpj|^(/i(P))). 

Since the Buchberger algorithm preserves the P[m] -homogeneity, we have 

Proposition 3.8 Let N be a left Ad+n[to]-submodule of Ad+n[toY generated by P[m]- 
homogeneous operators. Then there exists a Grobner basis with respect to -<h of N con- 
sisting of F[m]-homogeneous operators. Moreover, such a Grobner basis can be computed 
by the Buchberger algorithm. 



The following proposition can be easily proved in the same way as Theorem 3.12] 

Proposition 3.9 Let N be a left Ad+n-suhmodule of A^^j^^ generated by Pi, . . . ,Pi G ^d+„- 
Let us denote by h{N) the left Ad+n[to]-submodule ofAd+n[toY generated by h{Pi), . . . , h{Pi). 
Let G = {Qi(to), . . . ,Qfc(^o)} be a Grobner basis of h[N) with respect to -<h consisting 
of F[m]-homogeneous operators. Then G{1) := {Qi{l), . . . ,Qk{l)} is an F[m]-Grdbner 
basis of N. 

Thus we have an algorithm of computing an F[m]-Grobner basis for an arbitrary shift 
vector m G Z*". We can prove the following in the same way as p9l Proposition 3.11] 



Proposition 3.10 Let N andQjlto) be as in Proposition \3. 9{ andputM := Vx^Aj_+„N C 
T>x- Then for any germ P of M at p E Y , there exist germs Uj of T>x at p such that 
P = f/iQi(l) + ■ ■ ■ + UkQkil) and ordy[m](f/,g,(l)) < ordy[m](P) for j = 1, k. 

If the leading exponent of P G A(i+„[to]'^ is \expfj{P) = (A, a, i) G N x L x {1, . . . , r}, 
we define the leading position \pfj{P) of P by i. For a,P & L and i G {1, .■■,r}, we put 

(a, i) V {P, i) := (max{ai, max{a2d+2n, /32d+2n}, i), 

(a,z) + (^, i) := {a + (3,i). 

Let N and Qjito) be as in Proposition |3]^ and put A := {(«, j) | 1 < « < j < 
k, IpHiQM) = IpniQjiio))}- For {i,j) G A, let Sij{to), Sji{to) G Ad+n[to] be mono- 
mials such that 

lexpH{Sjiito)Qi{to)) = lexpHiSij{to)Qj{to)) = lexp^(Qi(to)) V lexp^((5j(to)), 
lcoeiH{Sji{to)Qi{to)) = Icoef h {Sij (to) Qj (to)). 

Then by the Buchberger algorithm, there exist F-homogeneous Uijiito) G Ad+„[to] so that 
we have 

k 

Sji{to)Qi(tQ) — Sij{tQ)Qj{tQ) = Uiji{to)Qi{to) 

1=1 

and either Uijiito) 7^ or else 

\e-x.Y>H{Uiji{to)Qi{tQ)) ~<H lexp^((5i(to)) V lexp((5j(to)) 
for each / = 1 k. 



The proof of the following proposition is similar to that of |29, Theorem 3.13]: 
Proposition 3.11 In the same notation as in Proposition \3.1(\ , the left A^^j^^-module 

Syz(gi(l), . . . , Qfc(l)) := {(f/i, f/fc) G A'd+n I f^iQi(l) + " " " + UMl) = 0} 
is generated by {Vij{l) \ {i,j) G A} with 

(i) (j) 

Vij{to) := (0, • ■ ■ , Sji{to), . . . , —Sij{to), . . . , 0) — (Uijiito), . . . , Uijk{to)). 



Now let us describe an algorithm for computing a free resolution of M which is adapted 
to the filtration Fy[m] (cf. Proposition p.3|) . Let iV be a left A^+n-submodule of 
such that M = A^^^^^/N. 

First, starting with a given m G Z^, let {Pi, ...,Prj} be an F[m]-Grobner basis of 
constructed as in Proposition p.9| . Put 

mi := (ordi7'[m](Pi), . . . , ordi.[m](P,,J). 



and define ipi : 



^d+n by 



IpliQl, ...,Qrr) ■= QlPl 



Qri Pr 



Then we get a set of generators of the kernel Keri/^i by using Proposition 3.11 



By the same procedure as above with A^, r, and m replaced by Ker?/^i, ri, and mi 
respectively, we obtain a homomorphism ip2 '■ ^d+n ~^ ^d+n ^^^^ lmip2 = Keiipi. In 
view of Propositions 3.101 and p. 11 , the sequence 

is exact for any A; G Z with m2 G Z^^ defined by 

m2 := (ordi;'[mi](V^2(l, • • • ,0)), . . . , ordF[mi](V'2(0, . . . , 1))). 

Proceeding in the same way, we can obtain a free resolution ( p.3|) which is adapted to the 
P[m] -filtration for any given / G N. 



4 The 6-function of a D-module 

Let Ai he a left coherent Px-niodule on X. We assume that a left y4(i+„-submodule N 
of is given explicitly so that M = T>x M holds with M := A^+„/iV. Set 

Af := Vx ®yid+n ^ ^x- We fix an arbitrary shift vector m = (mi, . . . , rrir) G Z^ and 
put 

gr^H(Ar) := P^[m](Ar)/P^i[m](Ar), 
gr^[m](M) := F^[m]{M)/ F^-^[m]{M). 

They are left gry (Px)-niodules. The P[m] -filtration and the associated graded module 
are defined also for y4rf+„- modules. Moreover we have 

gr^[m](A^) = g4{Vx)gr'y[m]iM). 

We put d := tidt^ -|- ■ ■ ■ -|- tddt^- This is the unique vector field modulo Fy^iVx) that 
operates on Jy I Jy as identity. Let 6' be a commutative variable corresponding to "d. 



Definition 4.1 The b-function (or the indicial polynomial) b{6,p) G K[6] of A4 along Y 
with respect to the filtration -Fy[m] at p G V is the monic polynomial b{6,p) G K[6] in s 
of the least degree, if any, that satisfies 

b{^,p)gT^[m]{M)p = 0. 

If such b{6,p) exists, Ai is called specializable along X at p. If A4 is not specializable at 
p, we put b{6,p) = 0. The global b-function b{9) of M. along Y is defined to be the least 
common multiple of b{6,p) with p running through Y. 

It is known that the specializability does not depend on the shift vector m while the 
6-function can depend on it (cf. [Q). It is also known that if Ai is holonomic, then 
is specializable at any p G X ([|T§,|g, [^). 

First, we reduce to the case r = 1. For each z = 1, . . . , r, let vTj be the projection of 
gry[m](T'^) to the i-th component and put 

gr°-[m](Ar)» := {P G gr°.[m](Ar) | 7rj{P) = for j = i + 1, . . . , r}. 

Note that gry[m](A/')^*Vs^y[j^](-^)'"* can be regarded as a left ideal of gry(T'x) by the 
projection to the i-th component. Then we get the following lemma: 

Lemma 4.2 Under the above notation, b{'d,p) is a generator of the ideal 

n {m n (grO.[m](Ar)«/gr?.[m](Ar)(-i)),) . 

i=l 

Let us now assume that the order -< satisfies 

(a, i) -< {a',j) if i < j for a,a' E L and i, j G {1, . . . , r}. (4.1) 

Let G be an F[m]-Grobner basis of with respect to ~<p defined by -< as in Section ^ 
Then I := gry[m](A/')^*Vg^Y[™](-^)*'*~^'' is generated by 

{7r,(ay[m](P)) I P G G, ay[m](P) G gr^ [m] (AT) W } . 

Our next task is to compute the intersection X fl Vyltidi, . . . , td^^]. For this purpose, 
we introduce commutative indeterminates v = {vi, . . . ,Vd) and w = {wi, . . . , Wd), and 
work with the ring y4rf+„[w,w]. For an element P of Ad+n of the form 

P= E , 

its multi-homogenization mh(P) G Arf+„[t>] is defined by 

mh(P) := J2 ■ ■ ■ v''d'-^'-''H^d'^x''d! 



with := min{z/j — \ a^^ap 7^ 0}. Let -<mh be an order on N"^ x N*^ x L 3 {p,cr,a) 
defined by 

(p, a, a) -<mh {p', cr', tt') if and only if |p + o"| < \p' + cr'l 

or else |p + o"| = \p' + (t'\, a < a' (4.2) 

with an arbitrary well-order < on L satisfying ( p.l| ) . Fixing an i G {1, . . . ,d}, we assign 
weight 1 to Wi, dt-, weight —1 to Vi, ti, and weight to all the other variables. An element 
of Ad+nivyw] is said to be multi-homogeneous if it is homogeneous with respect to the 
weight above for each i = 1, ...,d. Thus mh(P) is multi-homogeneous for any P G 
Put := S'lKi ■ ■ ■ Sd^,a for K = . . . , Hd) e Z'^ with Sij = if j > and Sij := 
otherwise. Let s = (si, . . . , s^) be commutative indeterminates. Assume that P G A^+n 
is multi-homogeneous. Then we have 

with some Q{si, . . . ,Sd,x, dx) G A„[si, . . . , Sd] and k G Z'^. We put 

i>{P){si,...,Sd) :=(5(si,...,Sd). 



Proposition 4.3 LetX be a left ideal of gTyiVx) ■ Let Gq be a finite subset of gTyiAd+n) 
which generatesT. Let Gi be a Grobner basis with respect to -<mh of the ideal o/^^+nb) 
generated by 

{mh(P) I P G Go} U {1 - ViWi I i = 1, . . . , d}. 

We may assume that Gi consists of multi-homogeneous elements since so does the input. 
Then the left ideal X fl Vyltidt^, . . . , tddt^] of VYltidt^, . . . , tddt^] is generated by 

6-2 := {^{P)it,dt„. . .,tddtj I P G Gi n Ad+n}. 

Proof: Let P be an element of Gi (lAd+n- Since P is multi-homogeneous and free of v, w, 
there exists k G Z*^ so that Q := ip{P){tidt^, . . . ,tddtj = S^P- By definition, P belongs 
to the ideal generated by mh(Go) and 1 — ViWi {i = 1, . . . ,d). Setting Vi = Wi = 1, we 
know that Q belongs to I. 

Conversely, let P be an arbitrary germ of XflPy [tidt^ , ■ ■ ■ , tddt^] at p G F. Multiplying 
P by a polynomial in x which does not vanish at p, we may assume P G An[tidt^ , ■ ■ ■ , tddt^]- 
In view of the definition of the multi-homogenization and the fact that mh(P) = P, there 
exists p G N*^ so that v^P belongs to the ideal generated by mh(Go). This implies that 
P belongs to the ideal generated by mh(Go) and {1 — ViWi \ i = 1, . . . ,d} since 

P = (1 - vPwP)P + v^w^P 

and (1 — f ^if^) belongs to the ideal generated by {1 —ViWi | i = 1, . . . , c?}. Set Gif\Ad+n = 
{Pi, . . . , Pfc}. Then by the definition of -<mh and Gi, there exist Qi, ■ ■ ■ ,Qk £ Ad+n so 
that 

P = QiPi + ■■■ + QkPk. 



Since Pi, . . . ,Pk are multi- homogeneous as well as P, we may assume that such is also 
the case with Qi, . . . , Qk- Hence there exist k,^^\ . . . , n^''^ E Zi'^ so that 



P = ij{P) = QiS,wij{Pi) + ■■■ + QkS,i.MPk)- 

This completes the proof. [| 

Let X be as in Proposition \i.3\ . Now we have obtained a set of generators G2 of 
X n X'y[ti9ti, . . . , tddtj]. We identify each tidt^ with Sj. Then from G2, we can compute a 
set of generators G3 of the ideal Xn(9y[s] of Oy[s] by eliminating dx by means of Grobner 
basis in the Weyl algebra (see e.g. ||2^ for details). Then it is easy to obtain a subset G4 
of K[x, 9] which generates the sheaf of ideals 

J -.= 1(1 Oy[0] = X n Oy[s] n Ovie] 

with 6 = Si + ■ ■ ■ + Sd again by Grobner basis in the polynomial ring. Let us denote by 
J the ideal of K[x, 9] generated by G4. Then b{9,p) is a generator of 

JpnK[9] = {Oy[9])j,J n K[9]. 

Our final task is to compute the 6-function at each point of Y by using the input (^4. 
This is achieved by primary decomposition. Let us state the method in a more general 
setting, where we replace the variable 9 by the variables s = (si, . . . , s^) for the sake of 
generality: So let J be an arbitrary ideal of K[x, s] whose generators are given. For each 
point p of Y, put 

BiJ,p) := iOY[s])pJ n K[s], 

which is an ideal of K[s]. Let J = Qi n ■ ■ ■ H Qi he a primary decomposition in K[x,s]. 
Then by the flatness of ((9y[s])p over K[x, s] we have 

B{J,p) = B{Qi,p)n---nB{Qi,p). 

Each ideal on the right hand side can be computed easily by the following: 

Lemma 4.4 Let Q be a primary ideal of K[x, s] and put 

Vy(Q) := {a; G r = fsT" I f{x) = for any f E Q H K[x]}. 

Then we have 

BiQ,p) -- 



QnK[s] zfpEYYiQ) 
K[s] tfpEY\YY{Q). 



Proof: First assume p ^ Vy((5). Then there exists a{x) E K[x\ such that a{p) 7^ 0. This 
implies that B{Q,p) = K[s]. Next assume p E Vy(Q) and b{s) E B{Q,p). Then there 
exists a{x) E K[x] so that a{x)b{s) E Q and a{p) 7^ 0. Suppose b{s) does not belong to 
Q. Then we have a{xy E Q with some j G N since Q is primary. This implies a{p) = 0, 
which is a contradiction. Thus we have B{Q,p) C Q fl K[s]. The converse inclusion is 
obvious. H 



Lemma 4.5 Let J be an ideal of K[x,s]. Then we have 

n B{J,p) = JnK[s]. 

Proof: Let J = Qi H ■ • ■ fl Q/ be a primary decomposition. Then by the preceding lemma, 
we have 

n B{j,p) = n l\B{Q,,p) 

p&Y peY j=i 

= f]{Qj n K[s] I Qj n K[x] ^ K[x]} 

= JnK[s] 

since Qj D K[s] = K[s] if and only if Qj n K[x] = K[x]. W 

Returning back to the ideal J of K[s, 6] generated by G4, we have only to apply 
Lemma with s replaced by 6. This gives us an algebraic stratification of Y so that 
b{s,p) is constant on each stratum as a function of p. Moreover, Lemma tells us that 
the global 6-function of Ad along Y is simply a generator of Jr\K[6]. Thus the algorithm 
is summarized as follows: 

Algorithm 4.6 (The 6-function of M := Vx M) 

Input: M = A^^j^^/N with an A^+^-submodule of A^^^, and m e 7/ . 

1. Compute a Grobner basis G of with respect to the order -<p that is defined 
through ( |3.2| ) by using m and an order -< satisfying ( |3.1| ) and ([4 .11). 

2. For i = 1 to r do 

(a) Gi := {7ri(ay[m](P)) \P eG, 7rj((Ty[m](P)) = for any j > i}. 

(b) Let Gil be a Grobner basis of the left ideal of Ad+n[v,w] generated by 

{mh(P) I P G G,} U {1 - v,w^ \ t=l,...,d} 



with respect to an order -<mh satisfying 

(c) G,2 ■■= {^(P) e Ar,[s] \PeGan 

(d) Compute Jj := {Gi2)nK[x, 9] first by eliminating d^, then eliminating S2, ■ ■ ■ , s'^ 
after substitution 9 = si + ■ ■ ■ + Sd and s'j = Sj for j = 2, . . . , d; here {Gi2) 
denotes the left ideal of An[s] generated by Gi2. 

3. The global 6-function b{9) of M is the generator of HUii^i ^ -^[^D- 

4. For i = 1 to r do 

(a) Compute a primary decomposition Jj = fljLi Qij in K[x, 9]. 

(b) Compute (generators of) Qij fl K[s] and Qij (1 K[x] for j = 1, . . . , /j by elimi- 
nation. 



5. For each p the local 6-function b{6,p) is the generator of the ideal 

r I- 

n n {Qij n K[e] I g{p) = for any g{x) G Qij H K[x]}. 
i=i j=i 

Let us remark on the coefficient field: Suppose that the input is defined over a subfield 
Ko of K. Then the steps 1-3 can be done over Kq instead of K and b{6,p) divides b{6) 
for any p However, the primary decomposition in the step 4 must be one in K[x, 6] 
not in Kq[x,6]. In fact, we need a primary decomposition over an intermediate field Ki 
with Kq d Ki <Z K so that h{9) factors into linear polynomials in -ft'i[6']. If, e.g., Kq is 
the rationals Q, such Ki is computable. Hence the primary decomposition in the step 4 
is certainly computable if the input is defined over Q in view of e.g., 0, |T^, and 
gives the local 6-function at any p G F = K^. 

As a special case where all the computation can be done over Kq, suppose that the 
ideal Jfl Kq[6] is generated by a polynomial which is a multiple of linear factors over Kq. 



Then the step 4 of Algorithm |4.6| can be computed over Kq and the step 5 is true for 
any p G K^; one can easily verify this by considering a projection of K to Kq. Note that 
this is exactly the case with the classical Bernstein-Sato polynomial (cf. |29|,||3^,[^ for 



algorithms) and Kq = Q by virtue of Kashiwara's theorem on the rationality [p^]. 

At this occasion, let us make a correction to [0]: Lemma 4.4 of |^ does not hold in 



general; we need field extension as explained above. This correction does not affect the 
rest of g. 

Example 4.7 Put X := K^ 3 (ti, ^2, y, z) and Y := {{ti, ^2, x,y, z) E X \ ti = t2 = 0}. 
Put M := A^/I with the left ideal / generated by 

ti-x^ + t2-y^ + + 3x^(9(1, dy - 2ydt^ + Sy'^dt^, - 2zdt^. 

Then the 6-function p) oi M. := T>x M along Y aXp eY = K^ is given as follows: 

h{s,p) = s(.s-^){s-\)(s- ^) (s + ^){s + I) 

x(^ + 4) (^ + 1) (^ + ^) (« + tI) (^ + 1) 

if p = (0, 0, 0) G Y] b{s,p) = s if p G {(x, y, z) | — = 0, y^ — z"^ = 0}; and b{s,p) = 1 
otherwise. 



5 Restriction of a D-module 

We retain the notation of the preceding sections. Put 

Vy^x ■■= Oy ®Ox T^x. 

Then T>y^x has a natural structure of (Py, T'x)-bimodule. Let M. = Vx ®Aa+n ^ be 
a coherent Px-module with a finitely generated A^+^-module M = A^^^^/N. Then the 
(D-module theoretic) restriction of to F is defined by 

My := Vy^x (^Vx ^ 



in the derived category of left I'x-niodules (see [T^ for the derived category and derived 
functors). 

In general, let £ be a module and Cj (j G Z) be additive subgroups of C such 
that tiCj C -Cj^i holds for i = 1, . . . ,d and j G Z. Then for any integer k, we define the 
Koszul complex associated with C, = {Cj}jf=z and ti, . . . , by 

^ A+d®z A Z'^ ^ A+d-i®z A Z^ ^ /:fc®z A Z'^ ^ 0, 

where S is defined by 

d 

(g) A ■ ■ ■ A CiJ = ^ tju ® e; A Ci, A ■ ■ ■ A Cj^. 

for a subset {ii, . . . , ij} of {1, . . . , c?} with the unit vectors ei, . . . , of Z*^. We denote 
this complex by IC'{C,[k],ti, . . . ,td). When Cj = C for each j, we also denote it simply 

j 

by /C*(£,ti, . . . ,td). Here we regard Ck+d-j® A 7/^ as being placed at the degree —j to 
be compatible with the cohomology theory. 

In particular, )C'(Vx,ti, . . . ^t^) is quasi-isomorphic to Vy^x in the derived category 
of right Pjc-modules. Hence we can identify My with the complex 

/C'(Px, ti, . . . , U) M = }C'iM,h,...,U). 

Our purpose below is to describe an algorithm to compute each cohomology group 
7i*(A^y) (for i = 0, —1, . . . , — (i since it is zero for other i) under the assumption that 
Ai is specializable along Y. Let b{6) be the global 6-function of Ai with respect to the 
filtration -Fy[m] with a given m G Z**. In what follows, we can replace b{9) by the local 
6-function b{6,p) in order to compute Aiy locally, i.e., on a Zariski neighborhood of p. 

Proposition 5.1 Let k be an integer such that b{k) ^ 0. Then the Koszul complex 
IC' {gTy[m]{A4)[k],ti, . . . ,td) associated with {gry[m](7Vl)}jgz is exact. 

We shall prove this proposition in a slightly more general situation. Let Ad := K[t]{dt) 
be the Weyl algebra on the variables t = {ti, . . . ,td) and define a filtration on it and the 
associated graded module by 

Fk{Ad) := { E a,,t^d'^ \ W - < k}, gr,(A,) := F,(A,)/F,_i(A,). 

Note that gT{Ad) := ®kezE,^k{^d) is isomorphic to Ad- In particular, we can regard 
gTQ^Ad) as a subring of Ad. 

Proposition 5.2 Let C = 0jgz A' graded gi (Ad) -module; i.e., assume gTj{Ad)Ci C 
for i,j G Z. Assume moreover that there exists a nonzero polynomial b{6) G K[6] 
which satisfies bi^ + i)Cj = for any j G Z with = tidt^ + ■ ■ ■ + tddt^. Let k be an 
integer such that b{k) ^ 0. Then /C(£,[/c], ti, . . . ,td) is exact. 



Proof: We argue by induction on d. First suppose d = 1. Then fC*{C,[k],ti) is the 
complex 

— > Ck+i — ^ Ck — ^ 0. 
Assume u G Ck+i satisfies tiu = 0. Then we have u = since 

= b{tidt, +k + l)u = b{dt,ti + k)u = b{k)u. 

On the other hand, there exists P E Ad so that b{tidti + k) = tiP + b{k). Hence for an 
arbitrary v G we conclude v G tiCk+i from b{tidt^ + k)v = 0. 

Now assume the proposition is true with d replaced by — 1. It is easy to see, as in 
the case of the usual Koszul complex (see e.g., [^, p. 188]), that }C'{C,[k],ti, . . . ,td) is 
quasi-isomorphic to the complex associated with the double complex 

/C'(£.[A; + l],ti,...,trf_i) 

i U (5.1) 
}C'{C.[k],ti,...,U-i). 

Let us denote by Cj and C" the kernel and the cokernel of : -^^j+i — ^ jCj. Then 
jC.' := ©jez '^'j and C" := ©j^z '^j are graded gr(y4rf_i) -modules. For any u G Cj, we have 

= b{tidt, + ■ ■ ■ + tddt, + J + l)u 

= bitidt^ H h td-idt^_^ + St^td + j)M 

= fe(ti9ti H h td-idta_, + j)u. 

On the other hand, for v E Cj, let be its residue class in C'j. Then we have 

= b{tidt, + ■■■ + Udt, + j)v = b{hdt, + ■■■ + U-idt,^, + ])v. 

Thus both C and C" satisfy the conditions of the proposition with d replaced by — 1. 
By the induction hypothesis, the complexes /C'(£^[fc]; ti, . . . ,td) and /C'(£^'[/c]; ti, . . . ,td) 
are exact. Hence the vertical chain map of (|5.1| ) is a quasi- isomorphism, which implies 
that )C*{C,[k],ti, . . . ,td) is exact. [] 

Under the assumption of Proposition p.l| , we have b{i) + j)gry [m](A^) = for any 
j G Z. In fact, for P G gr^('Dx), we easily get b{-t} + j)P = Pb{d). This yields 

b{d + j)^{[m\{M) = b{^ + j)gT{{Vx)gTY[m\{M) 

= gri.(Px)K^)g4[m](A<) = 0. 



Hence Proposition ^]1| is an immediate consequence of Proposition f).2\ 
Now for m G Z"^, we define the Fy [m]-filtration on Vy^x by 

F^[m](D^^;,) := F^[m](D^)/(tiF^+Hm](I?^) + --- + t,F^+i[m](D^)) 

r 

~ {P = ^ ^ ai,p{x)d^d^ei \ a^p{x) = if |z/| > /c - m^} 

= e e 

«=1 |i/|<fc— rrii 



Theorem 5.3 For an arbitrary m G Z*", let us take a free resolution ( \2.1\) of Definition 
2.1 of length I = d + 1, which is adapted to the Fylm]- filtration. Let mi, . . . , m^+i be 
the associated shift vectors. Take two integers < ki so that the global b-function b{6) 
of A4 satisfies b{k) ^ for any integer k with k < k^ or else k > ki. Then Aiy ^■^ 
quasi-isomorphic to the complex 



F^^lmom^ 







(5.2) 



with ro = r and uiq = m, where ipj is a homomorphism induced by tpj . In particular, we 
have My = if b{k) ^ for any k e Z. 



Proof: For any A; G Z, the complex 



^F^[mo](DX')-0 



is quasi-isomorphic to Fy[m](7W) in view of the exact sequence (p.l|) of Definition 
Hence we know that /C'(gry[m](A1)[/i;], ti, . . . ,td) is quasi-isomorphic to the complex as 
so dated with the double complex 







IC'{gT'y[moWnklt,,...,t,) 



/C-(gr^[m,+i](I?5*+^)[A;],ti,...,t,) 



On the other hand, we have a quasi- isomorphism 

X:-(gr^K](P5)[fc],ti,...,t,) ^ gr^K](P?^^). 
Hence K.' {gry[m]{Ai)[k],ti, . . . ,td) is quasi-isomorphic to the complex 



gr^[m,+i](P;?l^^ 



g?y^[moWy\^] 



0. 



(5.3) 



By virtue of Proposition the complex ( [5 .31 ) is exact if b{k) ^ 0. This implies the 
theorem since we have 







for sufficiently small k & Z. [] 

Note that F^^ [mi](I?p^^)/F^»'^[m,](r'p^;^) is a free Dy-module of rank 

ri 

'^^{u eZ'^ \ ko- niij < \iy\ < ki - niij}. 

i=i 

Hence Theorem gives us an algorithm to compute each cohomology group 7i*(A^y). 
In fact, we have only to compute the cohomology groups of the complex ( ^.21) as left An- 
modules with Vy^x replaced by The flatness of Vy over An assures us that the 

generators of the cohomology group over An also give the ones over Vy- The algorithm 
is summarized as follows: 



Algorithm 5.4 (The cohomology groups of the restriction of A4 to Y) 
Input: M = A^^j^^/N with an Arf+^-submodule of 
Output: H-\M'y) = Vy (A^//*) for < i < rf. 

1. Choose an arbitrary m G Z^; e.g., take m = (0, . . . , 0) by default. 



2. Compute the global 6-function b{6) of Ai along Y by the steps 1-3 of Algorithm ^ 
with M and m as input. 

3. If b{6) = 0, then Al is not globally specializable along Y; quit. 

4. Let ko and ki be the minimum and the maximum integral root of b{6) = 0. If there 
is no integral root, then we have ?^*(Aly) = for all i; quit. 

5. Compute a free resolution 

.rd+i i>d+l .ra JPd .ri J^l ^ro _V ^ p, 

^d+n ^ ^d+n ^ ■ ■ ■ ^ ^d+n ^ ^d+n ^ ^ ^ 

of M adapted to the Fy [m]-filtration (cf. Proposition |2.3| ) and the shift vectors 
mi, . . . , m^+i by using Proposition |3.11| successively, or by using Theorem p. 10 . 



6. Compute the induced complex 



[m,+i] {An [dtY-^^^ ) [mo] {An [dt_ 







'"0 



as a complex of finitely generated free left A^-modules, where A„[c}t] is identified 
with Ad+n/{tiAci+n H h tdAd+n)- Put ?/'o := 0. 

7. Via Grobner bases of modules over An-, compute the — i-th cohomology group 
Ker V'j/Ini'j/'j+i of the above complex in the form A^^/Ii with a left A^-module /j for 
2 = d. 



Note that in the step 5 of the above algorithm, only ^ipi, . . . , V'io+i needed if one 
wants to compute only the —i-th cohomology groups for i = 0, . . . , io- In particular, one 
does not need the free resolution to compute only the 0-th cohomology as will be the case 
with Algorithm |6.4| . 

As a direct application of the algorithm above, we obtain an algorithm to compute 
the cohomology groups with coefficients in the formal power series solutions of Ai: 

Ext^^^{M,K[[x]]) (z = 0,...,n) 

under the assumption that A4 is specializable along Y := {0}. In fact, we can easily 
verify that there exists an isomorphism (see e.g. |]22| , p. 428] for the case K = C) 

Ext\,^{M,K[[x]]) - Ext],iM-y,K) ^ H'iM'y). 

If i^' = C and Ai is Fuchsian along Y in the sense of |^2| (this condition holds if Ai 



is regular holonomic in the sense of ||T9|), then by virtue of the comparison theorem, we 
have also an isomorphism 

where C{x} denotes the ring of convergent power series in x. Note that Kashiwara's 
index theorem (|jl6|, p. 127]) gives the local index 

5:(-l)MimcExt^^(A<,C{x}) 

i>0 

at G X in terms of some topological quantity associated with the characteristic cycle of 
M. 

Example 5.5 Let us consider M := A4/I, where / is the left ideal of A4 (with K = C) 
generated by 

x^da + X4d4 - ai, Xidi + x^da - a2, X2d2 + x^d^ - a^, 0184-8283, 

where 01,02,03 G C are parameters. Put At := Vx M with X = and Yi : = 
{{xi,X2,X3,X4) E X \ xi = X2 = X3 = 0}. The global 6-function of A4 along Yi is 
(s — 02) (s + oi — 02 — 03). Hence the cohomology groups of the restriciton of Ai to Yi all 
vanish unless 02 or Oi — 02 — 03 is an integer. If dl — Qj2 — = 0, we have by Algorithm 



5.4 



VyJVy^Xi8i {i = 0), 

-UU W \ - J i1^Yj'Dy,X484)^ {i = -1), 

« l-^^yj - ] VyJVy^X,8, {l = -2), 

{i< -3). 

The 6-function of Ai along the point Yq := {(0,0,0,0)} is s — 02 — 03. Suppose Oi = 
^2 = fls = 0. Then the cohomology groups Wi^Aiy^) of the restriction of Ai to Yq are 
C, C^, C^, C, for i = 0, —1, —2, —3, —4 respectively. Since Ai is regular holonomic, this 
implies that Ext^^(A^, C{a;}) is C, C^ C^ C, for i = 0, 1, 2, 3, 4 respectively. 



Example 5.6 Put X := 3 {x, y) and M:=Vx {M/I) with / being the left ideal 

generated by 

xd^ - x{xd^ + ydy + a) {xd^ + bi) , ydy - y{xd^ + ydy + a) {ydy + 62) ■ 
Then by the computation of the restriction of M. to (0, 0), we get 



Sxt\,(M,C[[x, 



rc (^ = 0), 

C2 (z = l), 
[C (i = 2). 



for generic parameters a, bi, 62 (this means that we perform the computation over the coef- 
ficient field K := Q(a, 61, 62))- In particular, we have J2'i=o{~^y dime Sxt^j)^{M., C[[x, y]]) 
0. On the other hand, the characteristic cycle of M is 

3{^ = 77 = 0} + 4{x - 77 = 0} + 4{t/ = e = 0} + {x - y = e + - 0} + 7{x = y = 0} 

as a cycle in the cotangent bundle T*X — {{x,y,^,r])}. Thus, by Kashiwara's index 
theorem we have 

2 

^(-1)* dime Ext {M, C{x, y}) = 3 - (4 + 4 + 1) + 7 - 1. 

Hence Ai is not regular at (0,0). 

dimcSxt^^^{A4,C[[x,y]]) = 1 implies that the system A4 admits one dimensional 
space of formal power series solutions at the origin. In fact, the (divergent) formal series 

E ^"^7:f;i;"^^^^^ ^-y", {cu := c(c + 1) . • • (c + m - 1) 

m,n=0 [^)m[^)n 

spans the solution space. 



6 Tensor product and localization 

In this and subsequent sections, we denote by X the afiinc space First let us describe 
an algorithm to compute the tensor product and the torsion groups of two holonomic 
"Dx-modules A4i and M.2- We suppose that left A„-modules A^i and N2 are given so that 

Mr.= r{X,Mi) = A::/Ni (z = l,2). 

Let TTi, 7r2 : X x X — > X he the projections to the first and the second component 
respectively and put 

Then the exterior tensor product is defined by 

Mi®M2 := VxxX ®V'^^^ {t^i^Mi ®k 'K2^M2)- 



First, let us describe this exterior tensor product more concretely. Let ui, . . . ^Ur^ be 
the residue classes of the unit vectors ei , . . . , e^^ of A^^ , and , . . . , f the residue classes 
of the unit vectors e'j^, . . . , e^^ of • Then as a P^xx -module, A1i®A^2 is generated by 
Ui ®Vj {1 <i < ri, 1 < j < ^"2). Let us denote by (x, y) the coordinate system of X x X. 
For P = (Pi, . . . , P^J G and Q = {Qi, Q^J ^ we write 

p ® g := (p,(x, ajg,(y, 9,)),, G p^^;;^^. 

Let X' be the left "D^xx-submodule of {l^'xxxY^^^ generated by the set 

{P ®e'j, ei^Q\P eNi, Q e N2, 1 < 2 < n, 1 < j < ra} 
and put I := T'xxx ®x>^^^ 

Lemma 6.1 Under the above notation, there is an isomorphism 
Proof: Put 

Then it is easy to see that X' C /C'. Hence we have a commutative diagram 

\ i 
where $ is a i^'-bilinear map defined by 

i j 

which is well-defined by the definition of I'. In view of the universal property of the 
tensor product, we know that the vertical map above is an isomorphism. This completes 
the proof. |] 

Hence M.i(SiM.2 is computable with Ni and N2 being given. Put A := {{x,y) G 
X X X \ X = y} and identify A and X by the map tti. Then by [|14| , Proposition 4.7], 
which obviously applies to algebraic "D-modules as well, we have 

Ml ®Ox -^2 ^ V^^XxX ®Vxxx {Mi®M2) 
= (A^i®A<2)a- 

Suppose that Aii and A^2 are holonomic. Then it is easy to see that A^i(8>A^2 is a 
holonomic P^xx-module since its characteristic variety is contained in the Cartesian 
product of those of A^i and A^2- Hence A^i®A^2 is specializable along A and the 
following algorithm is correct: 



Algorithm 6.2 (The tensor product and torsion groups of two X'x-niodules) 

hiput: Holonomic systems A^j = Vx ®a„ Afj, where Mj = A^/Ni with an A„-submodule 

Ni of for i = l,2. 

Output: Tor'^''{Mi,M2) = Ckior k = 0,...,n. 

1. From sets Gi of generators of N^, compute 

G3 := {P ®e'j, Ci ® g I P G Gi, g G G2, 1 < ^ < n, 1 < j < ra}. 

2. Let G4 be the result of the substitution yi = Xi + ti {i = 1, . . . ,n) for each element 
of G3; let J4 be the left ideal of = K[t,x]{dt,dx) generated by G4 with t = 
{ti, . . . ,tn) 

3. Apply Algorithm E]^ with A2n/l4 as input and d = n to obtain 
for /c = n. 



Example 6.3 Put A := and 

M := Vx/Vxxd.,, M := Vx/Vxx. 

First, the exterior tensor product is given by 

M®M = Vx^x/{Vxxxxd, + Vx^xy) 

with {x,y) G A X A. Its global 6-function along the diagonal is s, and by restricting 
M.®N to the diagonal we get 

Torf^{M,M)=Vx/Vxx (^ = 0,1). 

In the same way, we get 

Tm°''(VxlVx{xd, + l),VxlVxx) = ('i = 0,l). 

Let be a holonomic Px-niodule and let / G K[x] be an arbitrary non-constant 
polynomial. Then we immediately obtain an algorithm to compute the localization 
M.[f~^] := Ox[f~^] ®Ox by combining this algorithm with that of computing 
given in since Ox[f~^] is holonomic Theorem 1.3]). Since Ox[f~^] is flat over 
Ox, the higher torsion groups vanish. 

Algorithm 6.4 (The localization M[f~^]) 

Input: A holonomic system Ai = T^x ®a„ (^n/^) ^"^^ ^ non-constant polynomial / G 
K[x]. 

Output: M[f-^]. 



1. Compute the global Bernstein-Sato polynomial of / as follows (cf. [^), where 
s is a single indeterminate: 

(a) Letting t be a single variable and let I be the left ideal of A^+i generated by 
t - f{x) and d^^ + {df/dxi)dt for i = 1, . . . , n. 



(b) Let b(6) be that in the step 3 of Algorithm |4.6| with An+i/I as input and d = 1. 
Put bf{s) := b{-s-l). 

2. Compute a set of generators of J/ := {P{s, x, dx) \ An[s\ \ Pf^ = 0} by fS^, Theorem 
19]. 

3. Let u be the minimum integer root of bf{s) = and put := {P(z/, a;, S^,) | 
Pis,x,dx) e Jf}. Then Ox[r'] ^ Vx {A^/Jf{^))- 

4. Compute Ox[f^^] ®Ox which is obtained as the output of Algorithm ^]2| with 
An/Jf{v) and A'^^/N as input and i = 0. 



7 Algebraic local cohomology groups 

Let fi, . . . , G K[x] be arbitrary polynomials and put Y := {x E X \ fi{x) = . . . = 
fd{x) = 0} with X := K"-. Let be a coherent Px-module. Our purpose is to compute 
the algebraic local cohomology groups 7i*y](A^) with support Y defined by Grothendieck 
as Px-modules. Recall that 7i*y](A^) is defined as the k-th derived functor of the functor 

riY]{M) ■.= limnomox{Ox/JY;M), 

where J7y is the defining ideal of Y and the inductive limit is taken as m tends to infinity. 
Note that if is a holonomic "Dx-module, then so is 7i|y](A^) (cf. fl^, Theorem 1.4]). 

Put X := K'^ X X and identify X with the hnear subvariety {0} x X of X. We set 

z ■.= {{t,x)ex\u = f,ix) (z = i,...,rf)}. 

Then B[z] '■= T^[z]i^x) isomorphic to V^/X, where T is the left ideal generated by 

df- 

tj-f^x) {j = l,...,d), dx^ + Y.^dt, (^ = l,...,r^). (7.1) 

Let TT : X — > X be the projection and put 

A' ■.= {(t,x,y)eXxX\x = y}. 



Then we can identify X with A' by tt. In the same way as |]3T1 , Lemma 6.3] we get 



Lemma 7.1 

with 

where pi and p2 are the projections of X x X to X and to X respectively. 

In fact, this lemma follows from the fact that A' is non-characteristic with respect to 
B[z]^-^- The proof of PT| , Theorem 6.1] yields the following: 



Proposition 7.2 For any coherent T>x-module Ai, we have an isomorphism 
as left Vx -module for any i > 0. 

Algorithm 7.3 (Algebraic local cohomology groups 7Y[y](Al)) 

Input: Polynomials /i, • • • , /d G K[x\ and a holonomic "Dx-module M. = Vx (A^/N) 
with an A„-submodule N of generated by G. 

Output: nlY]{M) = d-d for i = 0, . . . , d with Y := {x e K"" \ fi{x) = ... = fd{x) = 0}. 

1. Let / be the left 74rf+2n-submodule of A^^_^2n generated by 

Gi := {{tj- fj{x))ek\l<J <d, l<k<r} 

d Qf. 

U (d,^ + E \l<i<n, l<k<r} 

j=l (J-^i 

u {P{y,dy)\PeG}, 

where the unit vectors of Al,. 

2. Apply substitution yi = Xi + Zi for z = 1, . . . , n to Gi and let the result be G2- Let 
J2 be the sub module of ^rf+2„ generated by G2- 

3. Compute the 0-th cohomology of the restriction of A^^j^2nl J2 to {(t, x, z) | 2; = 0} in 
the form ®Ai+,, i^d+n/ ^^s) ^y Algorithm Here we can assume /cq = ^1 = 
skipping the steps 1-4 of Algorithm p.4|. 



4. Compute A := W{{V^ ®a,^^ {A^d+J h))\,x^x) for -rf < 2 < 0. 



Finally, assume that Y is non-singular and let 6 : F — > X be the embedding. Then for 
a coherent Py-module J\f, 

is a coherent Px-module with support in Y; here we put Vx^y '■= T^x / ij^xfi + ■ ■ ■ + 
l^xfd), which has a structure of (Px, ^^y)-bimodule. Moreover, the functor i+ gives 
an equivalence between the category of coherent Py-modules and that of coherent Vx- 
modules supported by Y ([0, Proposition 4.2]). In terms of this equivalence, we can 
compute the cohomology groups of the restriction of a holonomic Px-module to F by 
using Algorithm |7.3| and the following isomorphism: 



Proposition 7.4 ([|14|, Proposition 4.3]) 

L+n\My) =ni^fiM). 

Example 7.5 Put K := C^^ 3 (x, y, z) and 

M := Vx/{Vxd^ + Vxdy + Vx{z^d, + z)). 

Then the local cohomology groups of Ai with support Y := {(x, y, z) & X \ xz = yz = 0} 
are given by 

r Vx/iVxx + Vxy + Vx{z^d, + 2z + l)) (z = 2), 
n\y]{M) = \ (z = l), 

[ Vx/iVxdl + Vxixd, - 1) + Vxdy + Vxz) (z = 0). 

Moreover, 7i?y, (A1) = Vxu is also isomorphic to 



Vx/{Vxd, + Vxdy + Vxz) = Vxv 



by the correspondence u ^— d^v and v xu, where u and v denote the residue classes of 
1 G Vx in respective modules. 

The localization of by z is given by 

M[l/z\ = Vx/{Vxd, + Vxdy + Vxiz^d, + z + l)). 

Note that the methods of or [^] cannot be appied since z : ^ M. is not injective. 



8 Functors in the analytic category 



The functors studied in the preceding sections have analytic counterparts (cf. Q'lHl)- 
Throughout this section, we assume K = C and denote by OJ^ and PJf the sheaves 
on X of holomorphic functions and of rings of differential operators with holomorphic 
coefficients respectively. For a left Px-module M = V'x/J^, we put M.^^ := V^®x>x 
For k E 7i, put 

F^(p-) := {P e V^Iy I P{J^y e {J^y-^ for any j > k}, 



where Jy^ is the defining ideal of Y in . Then for m = (mi,...,mr) G TI ^ the 
Fy [m]-filtrations are defined by 



where Ui, . . . ,Ur are the residue classes of the unit vectors of (V^y. The graded modules 
are defined in the same way as in Section |. Put Vy^^x ■= ^03^" ^3f > which is a 
DJf )"t>imodule. Then the restriction of Ai^'^ to Y is defined by 



The 6-function of Ai^'^ along y at p G F is defined to be the generator of the ideal 

{b{e) G C[9] I 6(79)gr^[m](A^"'^) = 0}, 
where is defined as in Section ^. 

Lemma 8.1 Let JV be a coherent leftVx-submodule ofT>x- Define the Fy[m\-filtrations 
on M and on M'^'^ as in Section ^ with a shift vector m. Then we have 

grO.[m](Ar-) = gTUVt)®grliv,)gT'y{rn]{U) 
= V^'®T,^.gi^y[m]{M). 

Proof: We can prove the first equality by the same method (considering syzygies in 
the graded module) as Theorem 3.16] (cf. also Lemma 1.1.2]), where the case of 



r = d = 1 is treated; the argument applies to this case with trivial modifications. The 
second equality follows from 

gY''y{Vx)=Vy[t^dt,,...MdtX g4{V''^)=V^-[t^dt„...,tM. 



Proposition 8.2 The b-function b^^{9,p) of Jvl^^ and the b-function b{9,p) of Ai with 
the same shift vector m coincide for any p &Y . 



Proof: This follows from Lemma PTT] and the faithful fiatness of Oy^ over Oy (cf. 
Lemma 4.4]). [] 

Proposition 8.3 We have for any z G Z, 



Proof: Since we can regard 

we have an isomorphism V^_^x — ^y' T^y~>x as {T>^, X>x)-bimodules. Combining 
this with the flatness of PJf over Dx, and that of Vy^x over Vy-, we get 



(A^-)^ = A^) 



<an\« nan (T^s^o 

X 



®vy ipY^x M) 

®Vy M'y. 



This imphes the assertion since is flat over Vy- [] 
Proposition 8.4 For Vx -modules M. and M , we have 

Proof: The assertion follows from Proposition ^.3| and |]14|, Proposition 4.7]. [] 

For an algebraic set Y of X, the algebraic local cohomology groups of A4^^ are defined 
to be the derived functors of the functor 

TyY]{Mn := lim nomo^^m'/{JrT'-,Mn. 



Proposition 8.5 For a left T>x-module M., a polynomial f G C[x], and an algebraic set 
Y, we have 

M^^ir'] = v^^ (^v, Mir'], 

Proof: The first equality is an immediate consequence of Proposition ^.4| . The sec- 
ond equality follows from Proposition 8^, and Proposition 7^ together with its analytic 
counterpart. [] 



9 Homogenized Weyl algebra and Schreyer's method 
for free resolution 

In this section, we work in a framework more general than is needed in the preceding 
sections. Let An = K[x]{d) be the Weyl algebra over a field K of characteristic zero 
with X = {xi, . . . ,Xn), d = {di, . . . ,dn), di = d/dxi. We introduce a weight vector 



w = {wi, . . . , Wn] Wn+1, • • • , W2n) G Z^" \ {0} that Satisfies Wi + Wn+i > for i = 1, . . . , n. 
For each integer i/ G Z, we put 

n n 

■■={?= a^px'^d'^ G A„ I = if ^ WiUi + ^ if^+iA > '^j, 

o,/3eN" i=l i=l 

where aap G K and the sum with respect to a,/? G N" is finite. For a nonzero P G v4„, 
let ord^(P) be the minimum integer k such that P G It is easy to see that 

ord^(P(5) = ordu;(P) + ord^((5) holds for nonzero P,Q E An. More generally, for a shift 
vector m = (mi, . . . , mj.) G Z'', we define a filtration P^[m] of AJ^ by 

i=l 

where the canonical generators of A^- For a nonzero P G A^, we put 

ord^[m](P) := min{fc G Z | P G P^[m](A;)}. 

Now we introduce the homogenized Weyl algebra, which was introduced from the 
second version (1994) of kan/sml [ p8| : 

Definition 9.1 Let A^^^^ be the algebra over K generated hj h, x = (xi, . . . ,Xn), and 
d = {di, . . . , dn) which satisfy the relations 

j j i — ^ 5 ^ j ^3 — ^ ' ^ i ^3 ^3 ^ — ^ ^3 ^ ' 

hxi — Xih = 0, hdi — dih = (i, j = 1, . . . , n), 

where 5u = 1 and Sij = if i 7^ j. We call A^'^^ the homogenized Weyl algebra. The 
substitution h = 1 defines a i^"-algebra homomorphism 

p : Al^^BP^ P\h=i G A„. 

An element P of A^^^ is uniquely expressed as a finite sum 

AeN,Q,/36N" 

with aAo-/? G -ft'. The total degree of P G A^'^^ is defined by 

deg(P) := max{A + \a\ + \(3\ \ axap 7^ 0} 

if P 7^ and deg(P) = -00 if P = 0. 

Now let us take another vector n = {rii, . . . ,nr) G Z'', which describes the shift with 
respect to the total degree. For P = (Pi, . . . , P^) G {A^^}'', we put 

deg[n](P) := max(deg(Pj) +ni). 

l<i<r 

Definition 9.2 1. An element P = E[=i E(A,a,/3)GL axafSih^x^d^^ei of {A^^^Y is said to 
be /;,[n] -homogeneous if there exists G Z so that axa/si = unless A + |a| + +ni = 
k. 



2. For P = I]i=i J2{a,i3)€L o.ai3iX"d^ei G AJ^, we define the /;,[n]-homogenization /i[n](P) G 
(AW)'^ by 

i=l (a,f3)GL 

with /c := max{|a| + + rij | a^/ji 7^ 0}. This is /;,[n]-homogeneous. 

If n is the zero vector, we denote h[n]{P) simply by h{P). 

Lemma 9.3 For P G AJ^ and Q G An, we have p{h[n]{P)) = P and h[n]{QP) = 
hiQ)h[n]iP). 

Definition 9.4 Let ^ be a monomial order (i.e. an order satisfying ( |3.1D ) on Lx{l, . . . , r} 
with L := N^". We denote by lexp_,(P) the leading exponent of P G with respect 
to -<. Then -< is said to be adapted to the filtration Pu;[m] if {w,a) + rrii < {w,l3) + rrij 
implies -< j) for a,P E L and i,j G {1, . . . ,r}, and if lexp^(ej) -< lexp^^XjdjCi) 

for any 1 < i < r and I < j < n; here we write {w, a) = WiUi for a = (ai, . . . , a2n)- 

We fix a monomial order -<onLx{l,...,r} that is adapted to the Ptt,[m] -filtration. 
Then we define an order on N x L x {1, . . . , r} by 

(A,a,i) ^h[n] j) if and only if 
A + |a| + nj < /i + |/?| + Uj or else 
X+\a\+ni = fi+ \P\ + nj, {X,a,i) -< (/i, 

for A, A' G N, a, G L and i,j G {1, r}. Then it is easy to see that ^h[n] is a well-order. 
For a nonzero element 

^ = E E «Aa;3^/i^a;'^5^e, (9.1) 

i=l X,a,l3 

of (A^^'^))'', its leading exponent (Aq, Oq, /3o) "^o) = l6xp^[j^](P) G N x L x {1, . . . , r} is 
defined as the maximum element of {{X,a,P,i) \ axapi 7^ 0} with respect to -<h[n]- Then 
the leading position lp^jjj](P) and the leading coefficient lcoefh[n](P) are defined to be iq 
and aAo,ao,/3o,«o respectively. We denote them simply by lexp(P), lp(P), and lcoef(P) if 
there is no fear of confusion. The following lemmas follow easily from Definition and 
the definitions of A^^'^ and ~<h{a]- 

Lemma 9.5 For P e {A^n'^Y and Q G A'^^^ we have lexp(QP) = lexp((5efc) + lexp(P) 
with k = lp(P). 

Lemma 9.6 If P E (A^'^Y is h[n\-homogeneous and Q G A^J^^ is h[0]-homogeneous, then 
QP is h[n\-homogeneous. 



Lemma 9.7 Let w : N x L x {1, . . . , r} L x {l,...,r} be the projection. Then 
lexp^(p(P)) = -a7(lexp^j„](P)) holds if P is h[n]-homogeneous. 

In view of the above lemmas, we can define tlie notion of Grobner basis in the ho- 
mogenized Weyl algebra and can employ the Buchberger algorithm, which preserves the 
h [n] -homogeneity : 

Definition 9.8 Let be a left ^(f)-submodule of (A^f))''. Then a finite subset G of 
is called a Grobner basis of N with respect to -<h[n] if 

E{N) := {lexp(P) \ PeN\ {0}} = |J (lexp(P) + L). 

PeG 

Note that G generates N ii G is a Grobner basis of N since -<hin] is a well-order. 

Proposition 9.9 Let N be a left An-submodule of generated by Pi, .... Pk. Let 
h[n](N) be the left A[^^ -submodule of {A'^'>Y generated by /i[n](Pj) for i = 1, . . . , /c. Let 
G = {Qi, • • • , Qs} be a Grobner basis of h[n]{N) with respect to -<h[n\- Then for any P G 
N, there exist Uj G An such that P — Uip{Qi) + ■ ■ ■ + Usp{Qs) and ord^[m](L'jp((5j)) < 
ord^[m](P) forj = 

Proof: There exists z/ G N such that /i''/i[n](P) belongs to h[n]{N). By the divi- 
sion algorithm in (A!^^Y, we can find /i[n]-homogeneous Ui,...,Us G (AC'))'' such that 
/i^/i[n](P) = C/iQi + • • • + UsQs and leMUkQk) ^h[n] lexp{h''h[n]{P)) for A; = 1, . . . , s. 
Applying the ring homomorphism p, we get P = p{Ui)p{Qi) -|- • • ■ -|- p{Us)p{Qs) and 
lexp^(p(f/,,)p(Qfc)) = tJ7(lexp(f/feQ,,)) ^ w{\ex^{h''h[n]{P))) = lexp^(P). Since Ui and 
Qi are /i[n]-homogeneous, this implies OTdyu[m]{p{UiQi)) < ord^[m](P). [| 

We use the same notation as in the preceding proposition. Put A := {{i, j) \ 1 <i < 
j < s, Ip(-Pi) = Ip(-Pi)}- Then for (i, j) G A, let Sij, Sji G A!^^ be monomials such that 

lexp{SjiPi) — lexp{SijPj) — lexp(Pj) Vlexp(Pj), Icoef (^jiPj) = Icoef (5'jjPj). 

By the Buchberger algorithm, there exist /i[n] -homogeneous Uijk G A^^ so that we have 

s 

k=l 

and either Uijk 7^ or 

lexp{UijkPk) -<h[n] lexp(Pi) Vlexp(Pj) 

for each k = l,...,s. The following is an analogue of F.O. Schreyer's theorem for the 
syzygies in the polynomial ring: 



Theorem 9.10 In the notation above, let -<' be the order on N x L x {1, . . . , s} defined 
by 

{a, fi) -<' {P, v) if and only if lexp(P^) + a -<h[n] lexp(P;^) + (3 
or else lexp(P^) + a = lexp(Pj,) + [3 and /i > z/ 

for a, (3 E ISi X L and /i, z/ G {1, . . . , s}. Put 

m' := (ord^[m](Pi), . . . , ord^[m](P,)), n' = (deg[n](Pi), . . . , deg[n](P,)). 

Then -<' is a well-order adapted to the F.^[m']- filtration and called the Schreyer order 
induced by -<h[n]- 

1. For {i,j) G A; 

(i) U) 

Vij . (0, . . . , Sji^ • • ■ 1 ^iji • • • 5 0) ij^ijli ■ ■ ■ 1 Uijg) 

is h[n']-homogeneous and {Vij \ G A} zs a Grobner basis with respect to -<' of 
the module 

Syz(Pi, . . . , P,) := {(f/i, f/J G [Al^^r I f/iPi + ■ ■ ■ + f/,P, = 0}. 

2. Put 

Syz(p(Pi), . . . , p(P,)) := {{U,, Us) G I U,p{P,) + ■■■ + U,p{P,) = 0}. 

Then for any P G Syz(p(Pi), . . . , p{Ps)) H F^[m']{AY) with u eZ, there exist Qij G 
An such that P = 'E{i,j)eAQijPi^ij) (^"^d that QijpiVij) G P^[m'](A^). 

Proof: For a nonzero P G {A^^^'Y of the form ( |9.1|) , we define its initial term by 

in(P) := aA«^,/i^x"e^e, eK[h,x,iY 

with (A,a,/3,i) = lexp(P), where = (^i, . . . ,^n) are commutative indeterminates. Let 
Sij be the monomial in K[h,x,C,] obtained by substituting ^ for d in Sij. Then we have 
Sjiin(Pi) - %in(Pj) = for G A. Now suppose U = {Ui, . . . , f/^) G Syz(Pi, . . . , PJ 
and put {Xo,ao, j3o,io) := maxi<j<s lexp(f/jPj). Define Uj G K[h,x,^] by 

UjCig = in(f/,eij if lexp(t/,P,) = (Aq, ao, Po, io), 

and put Uj = otherwise. Then J2j=i Ujin^Pj) = holds. By the definition of -<' we have 

\exp^,{Vij) = {lexp{sji),i), lexp^,(?7) = lexp^,((Mi, . . . ,Us)). (9.2) 

On the other hand, since (0, . . . , Sji, 0, . . . , —Sij, 0, . . . , 0) for {i, j) G A constitute a Grobner 
basis with respect to -<' of the syzygies on in(Pi), . . . , in(P<j) by virtue of Schreyer's the- 
orem for the polynomial ring (cf. [^, Theorem 15.10]), we know that 

lexp^,((Mi,...,u,)) G U ((lexp(sji),i) + (N X L)). 



This completes the proof of the first assertion in view of ( |9.2| ). 

The second assertion follows from the first and Proposition ^]9| |] 

Let iV be a left A„-submodule of generated by Pi, . . . , P^. Let h{N) be the left A^^^- 

submodule of {Aj^^Y generated by h(Pi), . . . , /i(Pfe) (the homogenizations with n = 0). 

Starting with h{N) and n = m = 0, apply the first part of Theorem |9.10| repeatedely. 

Then we get an exact sequence 

...J^ ^Al^)y2 ^ (^W)n ^ ^j^{h)y, ^ (AW)7/,(iv) (9.3) 
with Tq := r. Put mo = Hq = and 

rwi := (ord^[mi_i](p(?/^j(l,0, . . . ,0))), . . . , ord^[mi_i](p(?/^i(0, . . . , 0, 1)))), 
n, := (deg[n,_i](V^,(l,0,...,0)),...,deg[n,_i](V^i(0,...,0,l))). 



Applying the homomorphism p to (|9.3| ), we get an exact sequence 

■■■'-^ Al- Al- Al- ^ Al/N — . 0. (9.4) 



Moreover, in view of the second part of Theorem 9.10, the sequence 



P^[mi](A-) P^[mo](A-) ^ Ftmj^^)KN n P^[0](AD) 



is exact for any k E Z. That is, the resolution (|9.4| ) is adapted to the P^-filtration. Fur- 
thermore, we can prove the following in the same way as its counterpart in the polynomial 
ring (IP, Corollary 15.11]) 

Theorem 9.11 By arranging the Grobner bases appropriately, we can construct the free 
resolution so that ip2n+2 = 0. 



It is known that there exists a shorter projective resolution. It is an open question to 
obtain shorter resolutions. 

Remark 9.12 1. If each element of the weight w is non-negative and the order -< 
adapted to P^[m] is a well-order, then the above construction can be done directly 
in A^ without the homogenized Weyl algebra. 

2. If the weight w satisfies Wi + Wn+i = for i = 1, . . . , n, then we can work in Anlxo]*" 



as in Section H instead of the homogenized Weyl algebra. Then Theorem 9.10 holds 



with the /i[n]-homogenization replaced by the P^„[m]-homogenization defined by 

r 

i=l a,/3GN" 

with k := min{(w, (a, /?)) + | a^f^i ^ 0} for P = ELi Eq,/3gN" a^pix'^d'^ei. 
Theorem p.ll| also holds in this case. 



In the computation of the free resolution described in this section, we can employ the 
method of La Scala and Stillman [2^, which computes the 'Schreyer frame' (initial terms 
of the resolution) first, then computes the resolution by a selection strategy or in parallel. 
We have implemented this algorithm in kan/sml. Most of the examples of this paper have 
been computed by using this implementation. 

Example 9.13 We give an example of the Shreyer resolution and explain the limits of 
our method (cf. Example 7.1 in |^) caused by the complexity. Let / be the left ideal in 



, X2, xs, 0:4, X5, xe, dx^-, dx2, dx3, dx^, dx^, dxg, ^1,^2, ts, dt^, dt2, dt^) 

generated by 

X4X2 - X5X1 + ti, X4X3 - XqXi + ^2, 3:5X3 - X6X2 + h 

x^dt^ + x^dt^ + dx^,-X4,dti + X6<9t3 + dx^, -x^dt^ - x^dt^ + dx^ 
-X2dti - X3(9t2 + dxi,xidt^ - X3(9i3 + (9^,, xi^i^ + X2<9t3 + dx^. 

We want to obtain the restriction of Ag/I along ti = ^2 = ^3 = 0, but we have not 
yet succeeded in getting the restrictions of all degrees except zero by using our algorithm 
and our implementation; we could get only a huge resolution. The hardest part in the 
computation of the restriction is to compute the quotient of the kernel over the image 
of the truncated complex (the step 7 of Algorithm |5.4| ). We cannot compute the higher 
order quotients by our algorithm and implementations. 

We homogenize the ideal with the variable s; 

X4X2S — X5X1S + ti, X4X3S — XqXiS + t2, X5X3S — XQX2S + ts 

x^sdt^ + X6sdt2 + dx^, -Xisdt^ + x^sdi^ + S^^, -x^sdt^ - x^sdt., + dx^ 
-X2sdt^ - X3sdt2 + dxi, xisdt^ - x^sdt:, + dxr^,xisdt2 + X2sdt.^ + dx^, 

which we denote by /'. We compute the Schreyer resolution with the weight matrix 





The reduced (non-/i-homogenized) Grobner basis of /' consists of 55 elements. The Betti 
numbers (i.e., ri,r2, . . .) of the Schreyer resolution are 630,3329, . . . where the reduced 
Grobner basis of syzygies by the Schreyer order among the 55-elements of the reduced 
Gronber basis of /' consists of 630 elements. The syzygy among 3329 elements by the 
Schreyer order could not be obtained, because of the memory exhaustion, on three MMX 
Pentium PCs (166 MHz) with 64 Mega bytes of memory under the Linux operating 
system (version 2.0.30) with kan/sml (version 2.980129) and the open sml module for 
communication between distributed processors The first Grobner basis of /' can be 
obtained in two seconds. 



The number of Grobner basis is smaller in this case if we compute in the homogenized 
Weyl algebra. We homogenize the ideal / in the homogenized Weyl algebra with the 
variable h; 

X4X2 - X5X1 + tih, X4XS - XqXi + t2h, 0:5X3 - XeX2 + hh, 

x-,dt^ + xe^tj + hd^^, -x^dt^ + x^dt^ + hd^^, -x^dt^ - x^d^ + hd^^, 

-X2<9ti - X3(9t2 + hd^^, xidt^ - xsdf, + hd^.^, Xidt^ + X2dt., + hd^^, 

which we denote by I^. We compute the Schreyer resolution with the weight matrix 

Xi 

( 1 -l\ 

li 1 1 )■ 

The reduced Grobner basis of consists of 44 elements which are less than those of 
V-homogenized Grobner basis. The basis can be obtained in 1.5 seconds by the same 
system. The Betti numbers of the Schreyer resolution are 506, 2422, . . .. 

We conjecture that in the homogenized Weyl algebra the minimal resolution exists 
and can be constructed by the algorithm in |]2D| since it is a graded algebra, but the 
boundary maps are not adapted to the filtration in general. For example, let us consider 
the following free resolution in the homogenized Weyl algebra 



where 



and 



— ^ {Arr 4'^ 4'V(pi,P2) 

aP = C{h,x,y,d^,dy), 

Vi((i,o))=pi, MiOA))=P2, M^) = {-P2,Pi) 



Pi = 2yhdr, + 3x'^dy,p2 = 2xd^ + 2,ydy + 6/i^. 



So, we cannot replace the Schreyer resolution by the minimal resolution. 

Finally, the authors note the iteration approach to get local cohomology groups as in 
may also improve the performance of our algorithms, but it is a future problem. 

During preparing this paper, we knew the paper Castro- Jimenez et al.. Homogenizing 
differential operators, 1997, preprint. The contents of this section may have an overlap 
with their paper. 
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